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In Ref. [l| we evaluated the radiative decay width cr//o(600) — * 77 under the assumption that a is described as 
a quark-antiquark state with flavor configuration Tin = (uu + dd)/y/2. For this purpose we utilized both local and 
non-local interaction Lagrangians describing the coupling of constituent quarks to the scalar field and which also 
allows a consistent, gauge-invariant inclusion of the electromagnetic interaction. Since a non-local model description 
of a quark-antiquark bound state is in our view unavoidable, we demonstrated in a pure quarkonium description that 
r<j^ 77 < 1 kcV for a mass of M a < 800 MeV. In Ref. Q we also stress that a further inclusion of meson loops will 
enhance the value of r ff ^ 77 , but an explicit calculation was not performed. Our main conclusion in Ref. [l[ was 
therefore that a dominant or pure quarkonia interpretation of the a does not allow for a full explanation of currently 
available data Q on r CT _ >77 . 

The preceding Comment by E. van Beveren et al. [|| claims that in [l[ we were mistaken on essentially three points 
on which we briefly elaborate in the following. 

Evaluation of the quark triangle diagram: In Ref. [l[ we stress that an accurate evaluation of the quark triangle 
diagram for a — ► 77 generates a term which in general causes destructive interference when compared to the corre- 
sponding 7T° — * 77 amplitude. This additional term vanishes under the peculiar condition M a = 2m q , where m q is 
the constituent quark mass. The original citation in Ref. [l| of the authors of preceding Comment with respect to 
this technical issue referred to this peculiarity. The analytical results both deduced in Ref. [1[ and Ref. [3| for r <T ^ 77 
now completely agree in the case of a local Lagrangian formulation. 

Nonlocal description of quark-antiquark bound states: To further illustrate the need for a nonlocal Lagrangian 
formulation we first refer to the Nambu Jona-Lasinio (NJL) model, which originally is given in local form (see e.g. 
Refs. [1, 0]). Regularization of loop integrals requires the introduction of a sharp cutoff A or a cutoff function. 
Independent of the precise form of the cutoff function the important point is that the cutoff A ~ 1 GeV has a well- 
defined physical meaning: it is related to the nonperturbative nature of the underlying and fundamental theory of 
quarks and gluons, QCD, and it sets the corresponding low-energy scale. Note however that the cutoff A, together 
with the precise form of the cutoff procedure is not included in the original NJL Lagrangian. Once a physical cutoff 
of the order A ~ 1 GeV has been introduced in an effective theory, it should be consistently included in all diagrams, 
including those that are ultraviolet (UV) convergent. A simple way to introduce this cutoff function already at the 
level of the starting Lagrangian is to render it nonlocal. This is explicitly done, for instance, in Ref. @, S, fic| - 
On a quantitative level NJL models with a proper introduction of the regularization procedure deliver an upper 
bound with r CT ^ 77 < 1 keV in a pure quarkonium interpretation [llj . This finding is consistent with the results of 
Ref. [l[ , although explicit numbers will depend on dynamical details and for example the explicit values for the a and 
constituent quark masses. 

Even on more general grounds, the QCD Bethe-Salpeter approach or QCD motived quark models based on bosoniza- 
tion of the QCD generating functional (for a review see [l2j and @,IH) show that a nonlocal interaction of a meson with 
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its constituents -the quarks- naturally emerges out of quark-gluon-dynamics. One might argue about the precise form 
of vertex functions and quark propagators, but the very fact that a nonlocal interaction arises seems undisputable. 

We therefore still argue that a non-local description of quark-antiquark bound states with a typical intrinsic scale of 
about 1 GcV will result in values for r CT ^ 77 below 1 keV, with explicit quantitative numbers depending on dynamical 
details and, trivially, on the mass of the a. Please note that most of the analyses now agree on a pole position of the 
a near (500 - i 250) MeV (see also note on scalar mesons 0). 

Meson-loops: Because of the large width of the a the coupling to rnr and KK followed by final state interaction 
will have a strong impact on the radiative decay width. In the Comment Q the authors deduce a net effect due to 
meson loops of about 40% of their total two-gamma width r CT ^ 77 w 3.5 keV. In Ref. [ll[ meson loops also contribute 
by about 50% but resulting only in r CT ^ 77 w 1.03 keV in total. A recent model dependent analysis of irir and 77 
scattering data [l3[ deduces a total 27 decay width of ~ (3.9 ± 0.6) keV, where the bulk part can be explained 

by rescattering. The direct or bare cr-pole contribution results in only r{J" „, w (0.13 ± 0.05) keV, actually in line 
with our results of Ref. [l| • Again, a model- independent estimate of meson loop contributions to the 27 decay width 
of the a seems presently not available. An analysis by Pennington [lij confirms the large value for the 77 decay width 
with r CT ^ 77 = (4.1 ± 0.3) keV, although Oiler et al. (HI or Bernabeu et al. [IB] deduce in their analyses smaller values 
of 1.8 ± 0.4 keV and 1.2 ± 0.4 keV, respectively. 

If a large 77 decay width of the /o(600) will be confirmed in future, our theoretical analysis shows that this 
result cannot be explained by the quark-loop contribution alone. Then we have two options: (i) Discard a dominant 
quark-antiquark interpretation of the /o(600), in agreement with many recent works [l7j; (ii) Argue that the meson 
loops generate the -by far- dominant contribution. In this case, however, it will be rather difficult to extract precise 
information about the nature of scalar states from 77 decays. 
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